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1 ABSTRACT 



2 The Japanese Space Agency's Hayabusa II mission is scheduled to rendezvous 

3 with and return a sample from the near-Earth asteroid (162173) 1999 JU3. 

4 Previous visible- wavelength spectra of this object show significant variability 

5 across multiple epochs which could be the result of a compositionally hetero- 

6 geneous surface. We present new visible and near-infrared spectra to demon- 

7 strate that thermally altered carbonaceous chondrites are plausible compo- 

8 sitional analogs, however this is a tentative association due to a lack of any 

9 prominent absorption features in our data. We have also conducted a series of 

10 high signal-to-noise visible- wavelength observations to investigate the reported 

11 surface heterogeneity. Our time series of visible spectra do not show evidence 

12 for variability at a precision level of a few percent. This result suggests two 

13 most likely possibilities. One, that the surface of 1999 JUS is homogenous and 

14 that unaccounted for systematic effects are causing spectral variation across 

15 epochs. Or two, that the surface of 1999 JUS is regionally heterogenous, in 

16 which case existing shape models suggest that any heterogeneity must be lim- 

17 ited to terrains smaller than approximately 5% of the total surface area. These 

18 new observations represent the last opportunity before both the launch and 

19 return of the Hayabusa 11 spacecraft to perform ground-based characteriza- 

20 tion of this asteroid. Ultimately, these predictions for composition and surface 

21 properties will be tested upon completion of the mission. 



22 Keywords: Asteroids; Asteroids, rotation; Asteroids, composition; Photome- 

23 try; Spectroscopy 
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24 1 Introduction 



25 In June of 2010 the Japanese Space Agency's Hayabusa mission successfully 



26 completed the first ever sample return from an asteroid (Fujiwara et al. , 2006 



S. Abe et al. , 2006 Nakamura et al. , 2011 , and references therein). The success 



of this mission has spawned a successor, Hayabusa II, which is scheduled to 
rendezvous with and return a sample from the asteroid (162173) 1999 JUS. 
This object was discovered by the Lincoln Near-Earth Asteroid Research (LIN- 
EAR) survey on 1999 May 10. It is an Apollo-class near-Earth asteroid with 
one of the lowest values of Af for any object with an absolute magnitude 
less than H=2Cp^ and is thus a highly favorable space mission target. The 
current Hayabusa II timeline proposes a launch in 2014 or 2015, rendezvous 
with the asteroid in 2018, and return of a sample to Earth in 2020. Like its 
predecessor this mission provides the opportunity to directly test the validity 
of ground-based analytical techniques with pre-encounter predictions for the 



38 physical and chemical properties of the asteroid (e.g. Binzel et al. , 2001b). 



Previous observations of 1999 JUS were conducted following its discovery in 
1999 and during its favorable apparition in 2007. Ground-based photometry 



ijietermined a light curve period of about 7.6S hours (Abe et al. 2008 Miiller 



et al. , 2011), which has since been confirmed with photometry from the 2012 



apparition (Kim et al. , 2012). A weighted average of space and ground-based 



observations in the mid-infrared suggest an effective diameter of 870 ± SO 



ipeters and an albedo of 6.954 ± 0.005% (Hasegawa et al. , 2008 Campins et 



al. 2009; Miiller et al. 2011). Prior to the present study, visible-wavelength 



spectra had been obtained in May 1999, July 2007, and September 2007 (Fig- 
ure [l| Table [T]). These data display marked variability in spectral slope and in 



http:/ / echo.jpl.nasa.gov/~lance/delta_v/delta_v.rendezvous.html 
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49 the presence of absorption features. The spectrum from 1999 shows UV ab- 

50 sorption short-ward of about 0.65 fim and a featureless neutral slope at longer 



51 wavelengths (Binzel et al. , 2001a). The July 2007 data, which suffer from the 

52 lowest signal-to-noise, show an overall red slope with a pronounced absorp- 



53 tion feature ~ 10% deep centered between 0.6 and 0.7 /im (Vilas, 2008). The 

54 September 2007 spectrum is a combination of data from two nights (Table 

55 [T| and is largely featureless with neutral slope (Vilas, 2008). Though these 

56 three spectra are broadly consistent with taxonomic classification in the C- 

57 complex, such pronounced variability, in this case at the level of ~ 10%, is 

58 highly unusual. 



59 This spectral variability has been attributed to a compositionally heterogenous 



surface (Vilas, 2008). The absorption feature between 0.6 and 0.7 /im in the 



July 2007 spectrum was linked to Fe^+ — )■ Fe^+ charge transfer transitions in 
iron-bearing phyllosilicates. The absence of this feature in the May 1999 and 
September 2007 spectra suggests a non-uniform distribution of these aqueous 
alteration byproducts. The possible existence of hydrated minerals on the 
surface of 1999 JUS is of considerable interest for the purposes of mission 
planning and the potential sampling of primitive Solar System material. Such 
surface heterogeneity could be inherited from the asteroid's parent body or 
could be exogenous in origin, as is believed to be the case for the "black 



69 boulder" seen by Hayabusa I on Itokawa (Fujiwara et al. , 2006) and for the 



70 dark material imaged by the Dawn spacecraft on 4 Vesta (Reddy et al. , 2012a) 



71 We present observations of 1999 JUS from the 2012 apparition, the last fa- 

72 vorable observing window before both the scheduled launch and return of 

73 the Hayabusa II spacecraft. Visible and near-infrared spectra were obtained 

74 to constrain the composition of the asteroid and to investigate claims of a 
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75 non-uniform surface. Other investigators obtained visible spectra during the 



76 2012 observing window (Sugita et aL, 2012 Vilas, 2012; Lazzaro et al. , 2013). 

77 These other data are not directly addressed here, though we note that none of 

78 these works found evidence for a 0.7 fim absorption feature or any repeatable 

79 spectral variability through multiple rotation periods. We also present here 

80 new visible-wavelength photometry to constrain the rotational phases of the 

81 asteroid accessed by the 2012 spectroscopic observations. Our observations 

82 and reduction techniques are presented in Section |2} Analysis of these data 

83 are presented in Section |3| We summarize and discuss the results in Section 

84 H 



85 2 Observations and Data Reduction 

86 Data presented here were taken with seven different instruments at six tele- 

87 scopic facilities and include visible and near-infrared spectroscopy, and visible- 

88 wavelength photometry. With the exception of one near-infrared spectrum 

89 from 2007, all data were obtained between October of 2011 and July of 2012. 

90 Tables [T] and |2] summarize the observing circumstances. 

91 2.1 Visible-Wavelength Spectroscopy: LDSS3 

92 New visible wavelength spectra were obtained with LDSS3 at the Magellan 

93 Clay 6.5m telescope on the nights of UT 2012 June 1-3, coinciding with the 

94 opposition of the asteroid at a phase angle of 0.2° on UT June 1 at 07:30. The 

95 instrument was operated with a 1.5" -wide and 4'-long slit and its VPH-AU 

96 grism. These settings produced a single-order spectrum with a useful spectral 

97 range of 0.44-0.94 fim at an average resolution of about 430 and a dispersion of 
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98 
99 
100 



1.89 A/pixel. Each individual spectrum represents a set of 4 x 180s exposures, 
split between two different nod positions approximately 20" apart along the 
slit. 



101 Solar analog stars were observed several times throughout each night for cali- 

102 bration purposes (Table [T|. With the exception of HD149182, the calibration 

103 stars were selected from a list of well-vetted analogs that were used extensively 



104 throughout the SMASS survey (Bus and Binzel, 2002). HD149182 was also in 



105 eluded because of its close proximity (< 10°) to the asteroid and was selected 

106 based on its photometric colors and G2V spectral classification. The high de- 



107 gree of consistency amongst our measured spectra (Section 3.2) confirms that 

108 HD149182 is a good solar analog. 

109 Reduction of these spectra employed standard IDL and IRAF routines, in- 

110 eluding a set of IRAF tools developed specifically for LDSS3. Combined He, 

111 Ne and Ar arc lamp spectra were obtained at the beginning and end of each 

112 night to provide dispersion solutions. Exposures of a white screen illuminated 

113 by a quartz lamp were used for flat field calibration. 

114 The extracted asteroid spectra were divided by each of the solar analogs (typi- 

115 cally ~ 10 individual observations of solar analogs throughout the night). The 

116 best solar analog for a given asteroid spectrum was chosen based on mini- 

117 mization of residual telluric features long- wards of ~ 0.7 /im. After division 

118 by an analog, individual points in the spectra were then rejected based on a 

119 0.25 boxcar sigma-clipping routine with a prescribed cutoff of 3 sigma. In 

120 all cases less than 10% of the data points were rejected and those that were 

121 removed were predominantly distributed at the extremes of the spectra out- 

122 side the useful spectral range of 0.44-0.94 fim. Generally, less than 10 pixels 

123 were removed from within this final trimmed spectral range. After cleaning. 
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124 the asteroid spectra were then re-binned by a factor of 13 pixels to produce 

125 a dispersion interval very close to 0.0025 /im per channel, similar to that of 



126 the 1999 spectrum (Binzel et al. , 2001a) and a factor approximately 3 times 



127 than that of the 2007 data (Vilas, 2008). The error bar on each binned dis- 



128 persion element was set to the standard deviation of the original 13 pixels. 

129 The data were finally normalized at 0.55 to produce relative refiectance 

130 spectra (Figure [T]). 

131 The wavelength-averaged signal-to-noise ratio (S /N) of these spectra vary from 

132 ~ 60 on the first night (A, B in Figure [T]), to ~ 30 on the second night (C), to 

133 ~ 20 on the third night (D, E, F). Since exposure times were held constant, 

134 this drop in S/N is likely caused by several effects. First, spectra A and B were 

135 taken at solar phase angles less then 1° (Table [T| and thus may have captured 

136 an opposition surge in brightness. Second, the moon was both waxing and 

137 increasingly close to the target throughout the run (Table [T]), thus causing 

138 a progressive increase in background flux. Third, the target's geocentric and 

139 heliocentric distances were both increasing, causing a drop in brightness across 

140 the three nights. 

141 2.2 Near-Infrared Spectroscopy: FIRE and SpeX 

142 Near-infrared spectra were obtained in June and July of 2012 with FIRE at 



143 ijthe Magellan Baade 6.5m telescope and in September 2007 with SpeX (Rayner 



etaL)|2003|) at NASA's Infrared Telescope Facility (IRTF; Table |T|). FIRE was 



145 operated in its high-throughput prism mode with a 0.8 x 50" slit. These set- 

146 tings produced single-order spectra at a resolution of approximately 400 from 

147 0.8 to 2.45 /im. Individual exposures were limited to 180s to avoid saturation 

148 of telluric emission features and of thermal emission from the instrument and 
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149 telescope at wavelengths longer than about 2.2 /im. Exposures were obtained 

150 in ABBA sequences with the object offset along the slit by 9" for each nod 

151 position. On UT 2012 June 5 a total of 4 x 180s exposures were obtained. On 

152 UT 2012 July 10, when the asteroid was nearly 2 magnitudes fainter, a total 

153 of 18 X 180s exposures were obtained. 

154 At the IRTF, SpeX was configured in its low resolution (R = 250) prism mode 

155 with a 0.8" slit for wavelength coverage from 0.8 to 2.5 fim. All observations 

156 were made with the telescope operating in a standard ABBA nod pattern 

157 with individual exposures of 120s. On UT 2007 September 18 a total of 40 

158 individual exposures were obtained, on UT 2007 September 20 a total of 50 

159 individual exposures were obtained. 

160 For both near-infrared spectrographs the slit mask was oriented along the par- 

161 allactic angle at the start of each observation to minimize the effects of atmo- 

162 spheric dispersion. As with the visible-wavelength spectra, solar analogs were 

163 observed to correct for telluric absorption and to remove the solar spectrum 

164 from the measured reflectance (Table [T]). Reduction of the SpeX data followed 



Sunshine et al. (2004), reduction of the FIRE data employed an IDL package 



166 [jiesigned for the instrument and based on the Spextool pipeline (Gushing et 



ah 2004). 



168 A composite visible/near- infrared spectrum of 1999 JUS is shown in Figure 



169 2 The asteroid is a C-type in the DeMeo et al. (2009) taxonomic system. 



170 The visible portion is a weighted average of the six LDSS3 spectra in Figure 

171 [ij The near-infrared portion is a weighted average of the SpeX and FIRE 

172 data, and thus represents a combination of data across multiple epochs. Other 

173 than some variability in slope, no significant differences were seen amongst 

174 the near-infrared spectra within the S/N of the data. 
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175 2.3 Visible-Wavelength Photometry: IMACS, Lulin, Tenagra II, and Bosque 

176 Alegre 

177 Photometric observations (Table [2]) were conducted with IMACS at the Mag- 

178 ellan Baade 6.5m telescope at Las Campanas Observatory in Chile, with a 

179 Princeton Instruments EEV Ik CCD at the Lulin Im telescope in Taiwan, 

180 with a SITe Ik CCD at the Tenagra II 0.81m telescope in Arizona, and with a 

181 Tektronics Ik CCD at the Bosque Alegre 1.5m telescope in Argentina. These 

182 data were taken on ten separate nights between October 2011 and June 2012. 

183 IMACS was operated using the //2 camera, which has a 27.5' field-of-view 

184 covered by a mosaic of eight 2k x 4k CCDs with plate scales of 0.2 "/pixel. 

185 We employed only one of the eight chips with a subraster about 400x400" in 

186 size. We used a Bessell R-band filter with exposure times of 60 seconds. These 

187 observations spanned approximately 3 hours per night over three consecutive 

188 nights. The Princeton Instruments CCD at Lulin Observatory has a roughly 

189 11' field of view with a plate scale of 0.516 "/pixel. These observations were 

190 conducted with a Johnson-Cousins R-band filter over the course of approxi- 

191 mately 2 hours on a single night. Exposure times were set to 300 seconds. The 

192 SITe Ik CDD at the Tenegra II telescope has a 14.8' field-of-view with a plate 

193 scale of 0.87 "/pixel. A Johnson- Cousins R-band filter was used. These ob- 

194 servations spanned approximately 4.5 hours per night over 4 non-consecutive 

195 nights with individual exposure times set to 360 seconds. The Tektronics Ik 

196 at the Bosque Alegre 1.5m has a plate scale of 0.33 "/pixel and a field-of-view 

197 of 5.7'. Observations were conducted with a Johnson R-band filter, spanning 

198 approximately 3.5 hours on two consecutive nights. Exposure times were 240 

199 seconds the first night and 200 seconds the second. 
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200 Reduction of these data employed standard IRAF procedures. Aperture pho- 

201 tometry was performed with aperture radii roughly 4 times the measured 

202 seeing and a background annulus equal to twice that separation. Some indi- 

203 vidual exposures were rejected based on contamination from nearby field stars. 

204 Photometric calibration was achieved for the IMACS data with observations 



205 of the RU 149 Landolt standard field (Landolt, 1992). All other photometry 

206 was performed differentially relative to at least one half dozen on-chip field 

207 stars and then scaled to match the calibrated IMACS data. A full table of the 
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collected photometry is included in the online Supplementary Data. 



209 3 Rotational and Compositional Analyses 



210 Our visible-wavelength photometry is used to construct a phase-folded rota- 

211 tional light curve of 1999 JUS. This provides the means to constrain the specific 

212 rotational phases accessed by the spectroscopic observations and to quantita- 

213 tively address the previously reported surface heterogeneity. The combined 

214 visible/near-IR spectrum is used to place preliminary constraints on compo- 

215 sition. 



216 3. 1 Rotational Light Curve 



Miiller et al. (2011) find a rotation period for 1999 JUS of 7.6S ± 0.01 hours 



218 and Kim et al. (2012) find a period of 7.625 ± O.OOS hours. We use these 



219 periods as a guide to phase-fold our broadband photometry (Figure [S]). The 

220 highest S/N data came from the three nights of observations at Magellan 

221 (Table |2|. The magnitudes of all other observations are scaled to the mean 

222 value of the photometrically calibrated Magellan data. As with previous light 
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223 curve measurements (Miiller et al. , 2011 Kim et al. , 2012), our data suggest a 



224 small peak-to-peak amplitude of about 0.2 magnitudes. Attempts have been 

225 made to exclude data points contaminated by background field stars, however 

226 the low galactic latitude of the asteroid (< 20° for all observations) certainly 

227 adds to the scatter in the measured photometry. 

228 The light curve in Figure|3]has been phase folded with a synodic period of 7.631 

229 hours relative to the first IMACS observation on UT 2012-04-05 at 06:40 (JD 

230 2456022.778). Light time corrections have been applied relative to this first 

231 IMACS observation. This light curve is wrapped such that the first and last few 

232 points are repeated at phases < and > 1. Robustly fitting a light curve period 

233 is difficult because the photometric errors are comparable to the amplitude of 

234 the light curve and our best quality photometry from IMACS does not span a 

235 full rotation period. A unique period is not found when employing the Fourier 



236 analysis of Polishook et al. (2012). Thus we attempt to manually constrain 



the period based on three features in the IMACS data: a subtle peak at phase 
= 0.2, a minimum at phase = 0.4 and another peak around phase = 0.6. The 
period is adjusted in increments of 0.001 hours until these three features are 
roughly aligned for all data. We find that these features are not aligned for 
periods within ±0.005 hours of 7.631 hours. This result is consistent with the 



242 range of periods suggested by Miiller et al. (2011) and Kim et al. (2012), but 



is not meant to be a more accurate determination. 



244 3.2 Non-detection of Spectroscopic Variability 

245 The phase folded light curve is used to provide rotational context for each 

246 of the new spectroscopic observations (Figure [3]). This phase folding demon- 

247 strates that rotational coverage was obtained over approximately 60% of the 
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248 asteroid's surface at a longitudinal sampling of ~ 45°. The spectroscopic ob- 



Binzel et al. 


2001a 


Vilas 


2008) 



250 included in Figure |3] because current best estimates for the rotation period 

251 (e.g. Kim et al. , 2012) still have an associated error that is large enough to 

252 preclude linking observations separated by several years. 



253 We compute reflectance ratios to investigate the possibility of a heterogenous 

254 surface (Figure |4]). These ratios are computed by median combining spectra A 

255 through F and then dividing that median into each of the individual spectra. 

256 Any deviations from a fiat line would indicate the presence of a heterogenous 

257 surface. Typical deviations are around 1%, with the largest (just under 4%) 

258 seen at the bluest wavelengths in spectrum C. However, all deviations fall 

259 within the S/N of the data and thus suggest that the regions on the surface 

260 accessed by our observations are homogenous to within a few percent. 



261 Previously reported spectroscopic variability around 0.8 fim (Lazzaro et al. 



2013) is not seen in our data set. This non-detection suggests the possibility 



263 raised by these authors, that it may be attributed to either instrumental opti- 

264 cal fringing or residual effects from the telluric water absorption band centered 

265 near 0.8 fim. 



266 The rotational context for these spectroscopic data provides a test of our 

267 observation and reduction techniques. Figure [3] shows that spectra A and F 

268 were obtained at nearly identical rotation phases. The fact that these spectra 

269 are indistinguishable (Figure |4]), but taken on different nights and calibrated 

270 with different solar analogs (Table [T]), is a nice validation of the consistency 

271 of our methodology. 
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272 3.3 Spectroscopic Surface Coverage 



273 Shape models of 1999 JUS based on photometry from 2007 have been inde- 



274 pendently developed by Kawakami et al. (2010) and Miiller et al. (2011). In 



275 both cases, the small amplitude of the light curve produces nearly spherical 

276 shapes. These models are based upon different data sets, but because of is- 

277 sues of low S/N, do not offer a single unique solution. Unsurprisingly then, 

278 neither shape model produces light curves that match our data well (Figure 

279 [3]). Clearly the development of a new shape model that fits photometry from 

280 2007 and 2012 is an important task for future work. Due to this ambiguity it 

281 is difficult to absolutely state the surface latitudes and longitudes accessed by 

282 the full ensemble of spectroscopic observations dating back to 1999. However 

283 we can compare model predictions for the spectroscopically accessed surface 

284 regions to illustrate the plausible spatial extent of any non-uniform surface 

285 features. 

286 The spectra of 1999 JUS from September 2007 are fully consistent with the 

287 new data presented here (Figure [T|. Any surface heterogeneity must then 

288 be restricted to regions accessed by the observations in May 1999, when the 



289 asteroid displayed UV absorption short-ward of about 0.65 /im (Binzel et al 



2001a) and in July 2007 when the asteroid showed a prominent absorption 



291 feature between 0.6 and 0.7 /xm (Vilas, 2008). These regions would have to 



292 be limited in extent to avoid violating the lack of absorption features in the 

293 other spectra. 

294 Table [Sjpresents the two sets of predicted sub-observer latitudes and longitudes 

295 for all visible-wavelength spectroscopic observations. Figure [5] illustrates the 

296 sub-observer points predicted by the two shape models. Though these sub- 
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297 observer points are not meant to be absolute and may change with future 

298 improvements to the shape model, their relative separations are informative. 

299 In particular, we are most interested in the distance that separates the sub- 

300 observer points for featureless spectra (September 2007, June 2012) from the 

301 sub-observer points for spectra that display absorption bands (May 1999, July 

302 2007; B and V in Figure |5]). We can roughly calculate the physical distance 

303 between sub-observer points by assuming a spherical body with a diameter of 



304 870m, a reasonable approximation for this asteroid (Kawakami et al. , 2010 



Miiller et al. , 2011 ). The physical distance between sub-observer points is then 



306 just the distance along the great circle that connects these points. 



307 According to the Kawakami et al. (2010) model, the observations from May 



308 1999 were taken at a sub-observer point about 180m away from the next closest 

309 point which corresponds to spectrum F (June 2012). And the observations 

310 from July 2007 correspond to a sub-observer point about 160m away from 



311 the September 2007 spectra. According to the Miiller et al. (2011) model, 

312 the observations from May 1999 were taken about 140m from the September 

313 2007 points, and the observations from July 2007 correspond to a sub-observer 

314 point about 200m away from spectrum D (June 2012). Taken as a whole, and 

315 independent of which shape model is employed, these calculations suggest 

316 that the reported spectroscopic variability would be a consequence of regional 

317 surface heterogeneity on spatial scales smaller than about 400m in extent 

318 (200m in radius). Regions of this size would occupy about 5% of the asteroid's 

319 total surface area. 

320 The ambiguity in shape models emphasizes that rotational information is not 

321 correlated over multiple years. This is indicated by the horizontal bands in Fig- 

322 ure|5j The sub-observer points for those spectra with absorption features could 
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323 lie anywhere along these bands. This provides some freedom for the spectrally 

324 anomalous regions to "hide" in the gaps between those sub-observer points cor- 

325 responding to featureless spectra. However, the sub-observer points for both 



326 the Binzel et al. (2001a) and Vilas (2008) spectra would have to be shifted by 



327 very specific amounts to reconcile the full suite of visible spectra. This issue 

328 is further complicated upon consideration of the sub-observer hemispherical 

329 footprints of our spectroscopic observations (see online Supplementary Figure 

330 A). These footprints indicate that we accessed a significant majority of the 

331 surface, leaving very little surface area in which spectroscopically anomalous 

332 regions could exist. 

333 3.4 Compositional Analogs 



334 With the full visible/near- infrared spectrum (Figure |2]) we can make prelim- 

335 inary statements about the composition of 1999 JUS. However, these state- 

336 ments must be interpreted with some skepticism due to the lack of promi- 

337 nent absorption features and issues of degenerate compositional interpreta- 



tions when dealing with such featureless spectra (e.g. Vernazza et al. , 2009 



Reddy et al. , 2012b|[c ). A preliminary attempt at compositional analysis is 



justified by the potential for predictions to be directly tested with a returned 



341 sample. We base this analysis solely on our new observations; see Vilas (2008) 

342 for an interpretation of previous visible-wavelength spectra. 

343 Our visible and near-infrared spectra confirm a C-type taxonomic classifica- 

344 tion, but show no evidence for the previously detected 0.7 micron phyllosilicate 



345 band (Vilas, 2008) or a prominent UV drop-off short-wards of about 0.65 /im 



346 (Binzel et al. , 2001a). The data may indicate weak bands around 0.9 fim and 



347 short-ward of 0.7 fim, but these features are near the limit of the S/N and 
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348 are not robust detections. Without any strong features to directly analyze 

349 we employ spectral matching techniques to constrain plausible compositional 

350 analogs. Comparison spectra were obtained from a July 2012 release of Brown 



351 University's RELAB database (Pieters, 1983) containing reflectance spectra 

352 for just over 17,000 samples. Our search began by selecting only those samples 

353 with a spectral range overlapping that of the asteroid and only those samples 

354 with an albedo within 5% of that measured for 1999 JUS. This narrowed the 

355 list of comparison spectra to 8024 samples. For each of these, either the as- 

356 teroid or the meteorite data were interpolated to match the wavelengths of 

357 the spectrum with the smallest dispersion interval. Both spectra were then 

358 normalized at the minimum asteroid wavelength (0.44 /zm) before a reduced 

359 difference was calculated. 



360 The two meteorite spectra with the lowest values are shown in Figure 

361 [2j These meteorites are a thermally altered sample of the CM carbonaceous 

362 chondrite Murchison and the unusual CI chondrite Yamato 86029. The reduced 

363 values for these meteorites were 2.89 and 3.25 respectively. These were 

364 clearly the best fits as the next lowest was 6.75 for a sample of a mid-ocean 

365 ridge basalt (RELAB ID RB-CMP-023) that was ground to particle sizes 125- 

366 250 ^m. This spectrum of oceanic basalt matched the general slope of 1999 

367 JUS but showed a subtle 1 /im absorption band and thus is unlikely to be a 

368 good analog. 



369 The mean albedos of the two best-fit samples (6% for Murchison, S% for Yam- 

370 ato) are by design of the search algorithm close to the ~7% measured for 1999 



371 JUS dHasegawa et al.[|2008[ [Campins et al.||2009t|Muller et al.||201lD . The two 

372 best-fit spectra both correspond to the smallest available particle size sorting, 

373 < 6S /im for Murchison and < 125 /im for Yamato. The Murchison spectrum 
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represents one in series of heating experiments where the temperature of the 



375 sample was raised in increments of 100°C from 300-1000°C (Hiroi et al. , 1996). 



At each temperature increment a reflectance spectrum was measured. Heat- 
ing was shown to reduce spectral slope and decrease the depth of absorption 
bands. The Murchison spectrum corresponding to a temperature of 900° C was 
the closest match to 1999 JUS. Yamato 86029 has unusual textural properties 
and a chemical composition indicative of a rare class of CI chondrites that ex- 



iperienced moderate thermal alteration at temperatures of 500-600° C (Tonui 
et al. , 2003[ ). 1999 JUS has been previously linked to carbonaceous chondrites 



383 (Vilas, 2008 Lazzaro et al. , 201S), the analysis here expands upon these re- 



sults by specifically identifying thermally altered samples as the best available 
spectral analogs. 



4 Discussion 



We have presented new observations of the Hayabusa II target asteroid 1999 
JUS taken during the favorable 2012 apparition, the last opportunity before 
both the launch and return of the spacecraft to perform ground-based char- 
acterization. Visible wavelength photometry was used to produce a rotational 
light curve. A period of 7.6S1 ± 0.005 hours and a small amplitude of a few 
tenths of a magnitude are consistent with independent light curve measure- 



393 ments (Miiller et al. , 2011 ; Kim et al. , 2012 ). Visible and near-infrared spectra 



394 were used to constrain taxonomy and composition. The asteroid is a C-type in 



395 the DeMeo et al. (2009) taxonomic system and the best spectral analogs are 



thermally altered carbonaceous chondrites. However, the lack of any promi- 
nent absorption features makes this compositional association tentative. High 
quality spectra that extend into the near-UV (O.S-0.4 /im) or the S /im region 
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399 could provide further compositional clues if diagnostic absorption bands are 



400 detected (e.g. Hiroi et al. , 1993 Rivkin et al. , 2002 ) and may serve as a test 



401 ipf the coUisional and thermal evolution of primitive asteroid types (Vilas and 



Sykes, 1996) 



The putative link to carbonaceous chondrites and particularly those that are 
thermally altered is interesting and merits some preliminary speculation. CM 
and CI chondrites are primitive meteorites with appreciable quantities of or- 
ganics and water of hydration (Cloutis et al. , 2011a|[b ). The possible return of 
such samples would be relevant to issues of astrobiology and to understanding 



primordial chemistry in the solar nebula (e.g. Morbidelli et al. 2000). The 
suggested link to heated carbonaceous chondrites could imply that surface 
temperatures on 1999 JUS were significantly higher at some point in its past. 
There is roughly a 50% probability that the chaotic orbital evolution of 1999 
JUS resulted in a perihelion inside of 0.5 AU (its current perihelion distance 
in 0.96 AU), which translates to surface temperatures of approximately 500 



K (Marchi et al., 2009 Michel and Delbo 2010). This is comparable to the 



in-situ alteration temperatures experienced by Yamato 86029, but is less than 
the 900°C experienced by the Murchison sample. Nevertheless, this highlights 
that thermally altered primitive meteorites could be relevant to understanding 
the composition of 1999 JUS. Ultimately, the results from Hayabusa II will 
provide the first-ever direct test of such compositional interpretations from 
largely featureless asteroid spectra. 



421 We have also presented a time series of visible spectra to address previously 



reported spectroscopic variability (Vilas 2008 Lazzaro et al. 201S ). Our spec- 



423 tra show no variability at the level of a few percent. There are several possible 

424 explanations for this non-detection. First, the asteroid's surface may be spec- 
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troscopically and compositionally uniform, in which case some unaccounted 
for systematic effects would be the cause of spectral variation across epochs. 
The hemispherical footprints of our spectral observations cover nearly the en- 
tire surface of the asteroid (Supplementary Figure A), thus supporting this 
hypothesis. Variations in spectral slope are common with asteroid observa- 
tions and can be attributed to observational issues such as weather, phase 



431 angle (Sanchez et al. , 2012; Reddy et al. , 2012d), or the use of imperfect solar 



analogs, but it is unclear how these issues would result in variable absorption 



(features. Other than a moderately high airmass of ~ 1.9 for the |Binzel et 
al. ( 2001a[ ) spectrum, we find no obvious observational reasons why the two 



435 spectra with absorption features should be different. 

436 Multi-epoch spectral variability could alternatively be linked to recent alter- 

437 ation events that induced global or near global changes in surface properties. 



438 Planetary encounters (Binzel et al. 2010) or impacts (Richardson et al. 2004) 



439 can alter reflectance properties, however these seem unlikely because two such 

440 events would be required in just 13 years to explain the observed spectroscopic 

441 variabihty from 1999 to 2007 to 2012. 



442 A final explanation is that the surface of 1999 JUS is in fact heterogenous. 

443 In this case spectral variability could result from the obscuration of surface 

444 regions by local topography as a consequence of changing viewing aspect. 

445 However, the orbital longitudes of both the Earth and the asteroid were very 

446 similar during the May 1999 and June 2012 observations. This would then 

447 require a highly irregular shape (i.e. large topographic relief) to result in sig- 

448 nificant obscuration and spectral variability. A highly irregular shape is in- 



consistent with current shape models (Kawakami et al. , 2010 Miiller et al 



450 2011). Therefore, spectral variability caused by a non-uniform surface would 
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have to be attributed to rotational effects. The combination of rotational data 
and time resolved spectra provides a means to address this possibility. Exist- 



453 ing shape models (Kawakami et al. , 2010 Miiller et al. , 2011) suggest that 



the maximum physical separation of sub-observer points for the full ensemble 
of spectroscopic data is about 200 meters (assuming that the asteroid is 870 
meters in diameter; Figure |5]). Therefore any surface heterogeneity must be 
confined to regions no larger than about 400 meters in extent, which corre- 
sponds to a size independent surface area fraction of around 5%. However, we 
note again that current shape models do not predict light curves consistent 
with photometry across multiple epochs. Therefore, the analysis here con- 
strains the plausible spatial extent of non-uniform surface features, improved 
shape models will be required to refine the spectroscopically accessed terrains 
indicated in Figure [5] 



464 Few asteroids, particularly sub-km objects like 1999 JUS, have been so ex- 

465 tensively observed across multiple epochs. However, non-uniform surfaces are 

466 generally believed to be rare. The long-time exception is asteroid 4 Vesta, one 

467 of very few objects where variability over multiple rotation periods has been 

468 confirmed with remote observations. Vesta displays regional and hemispheri- 

469 cal color heterogeneity, in some cases well in excess of several lO's of percent 



470 (e.g. |Bobrovnikoff[ |1929| [Gaffeyj |1997| [Reddy et aL| |2012a[ ). However, Vesta 

471 is much larger than 1999 JUS, is a different spectral class, and thus may not 

472 be the best object for comparison. 



473 There have been reports of rotational spectroscopic heterogeneity amongst 

474 C-class asteroids. For instance, 10 Hygiea, 105 Artemis, 1S5 Hertha, and 776 

475 Berbericia each show variability at levels less than about 5% in the depths of 



476 their 0.7 /im features (Rivkin et al. , 2002 Vilas, 2008). Detecting such subtle 
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477 variability is intrinsically difficult and as far as we are aware remains to be veri- 

478 fied for any asteroid over multiple rotation periods. Asteroid 21 Lutetia, which 

479 was visited by the Rosetta spacecraft, shows variability in visible- wavelength 



480 spectral slope at the level of a few percent (Barucci et al. , 2012). The S/N of 

481 our data for 1999 JUS may be insufficient to detect such muted variability. 

482 But again, these examples of spectral variability occur on asteroids that are 

483 much larger than 1999 JU3. 



The former target of the European Space Agency's Marco Polo mission (175706) 
1996 FG3 and asteroid 25143 Itokawa, the target of the Hayabusa I mission, 
are examples of well-studied asteroids in a size regime closer to that of 1999 
JU3. 1996 FG3 shows prominent (~ lO's of percent) spectral variability in the 



488 near-IR (de Leon et al. , 2011 ). The origin of this variability remains a mystery. 



489 Itokawa displays up to 15% color variability on its surface (Saito et al. , 2006 



M. Abe et al. , 2006). This non-uniformity is limited to regions on the order of 



lO's to 100 meters in extent or at most a few percent of Itokawa's surface area. 
This extent is similar to the constraints we have placed on 1999 JU3's surface, 
though it remains unclear whether such small regions could significantly in- 
fluence ground-based observations that only capture hemispherical averages. 
Ground-based observations of Itokawa prior to arrival of Hayabusa I did not 



496 exhibit any variability (Binzel et al. 2001b Abell et al. , 2007). 



497 For now, the issue of heterogeneity on the surface of 1999 JU3 remains an 

498 open question that may find some resolution from updated shape models. The 

499 arrival of the Hayabusa II spacecraft will ultimately determine whether the full 

500 ensemble of visible spectra are truly indicative of a non-uniform surface. More 

501 generally, the possibility of multiple planetary missions to near-Earth asteroids 

502 in the coming decades will provide an opportunity to further investigate the 



23 



503 extent of surface heterogeneity on sub-km asteroidal bodies, and thus provide 

504 clues about formational and evolutionary processes in the Solar System. 
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UT Date Instrument and Facility 



V A (AU) R (AU) a (°) 



2011- 


-10- 


-30 


EEV Ik CCD Lulin Im 


20.1 


0.438 


1 


.002 


76.1 


2012- 


-04- 


-05 


IMACS Magellan 6.5m 


19.8 


0.461 


1 


.237 


49.4 


2012- 


-04- 


-06 


IMACS Magellan 6.5m 


19.8 


0.457 


1 


.240 


48.9 


2012- 


-04- 


-07 


IMACS Magellan 6.5m 


19.8 


0.454 


1 


.243 


48.5 


2012- 


-05- 


-31 


SITelk Tenagra 0.8m 


17.7 


0.352 


1 


.366 


0.7 


2012- 


-06- 


-07 


SITelk Tenagra 0.8m 


18.2 


0.366 


1 


.377 


7.1 


2012- 


-06- 


-08 


SITelk Tenagra 0.8m 


18.3 


0.368 


1 


.378 


8.2 


2012- 


-06- 


-09 


SITelk Tenagra 0.8m 


18.4 


0.371 


1 


.380 


9.2 


2012- 


-06- 


-20 


Teklk Bosque Alegre 1.5m 


19.0 


0.409 


1 


.393 


19.4 


2012- 


-06- 


-21 


Teklk Bosque Alegre 1.5m 


19.0 


0.413 


1 


.394 


20.3 



Table 2: Summary of 1999 JU3 photometric observations. 
The columns in this table are the same as in Table [U 
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UT Midpoint Reference V A (AU) R (AU) a (°) Sub-Observer Long, and Lat. 



Miiller et aL (2011) Shape Model 



J L 



1999-05-17 07:49 
2007-07-11 10:05 
2007-09-10 07:00 
2007-09-11 05:50 
2012-06-01 07:07 
2012-06-01 08:09 
2012-06-02 03:54 
2012-06-03 01:54 
2012-06-03 03:34 
2012-06-03 04:54 



BOl 
V08 
V08 
V08 

A 

B 

C 

D 

E 

F 



17.7 
20.3 
17.9 
17.9 
17.7 
17.7 
17.8 
17.9 
17.9 
17.9 



0.305 
0.569 
0.273 
0.271 
0.353 
0.353 
0.354 
0.356 
0.356 
0.356 



Kawakami et al. 



(2010) Shape Model 



1.314 
1.381 
1.254 
1.251 
1.367 
1.367 
1.368 
1.370 
1.370 
1.370 



6.1 
40.3 
22.5 
22.8 
0.2 
0.2 
1.0 
2.0 
2.1 
2.2 



1999-05-17 


07:49 


BOl 


17.7 


0.305 


1 


314 


6.1 


2007-07-11 


10:05 


V08 


20.3 


0.569 


1 


381 


40.3 


2007-09-10 


07:00 


V08 


17.9 


0.273 


1 


254 


22.5 


2007-09-11 


05:50 


V08 


17.9 


0.271 


1 


251 


22.8 


2012-06-01 


07:07 


A 


17.7 


0.353 


1 


367 


0.2 


2012-06-01 


08:09 


B 


17.7 


0.353 


1 


367 


0.2 


2012-06-02 


03:54 


C 


17.8 


0.354 


1 


368 


1.0 


2012-06-03 


01:54 


D 


17.9 


0.356 


1 


370 


2.0 


2012-06-03 


03:34 


E 


17.9 


0.356 


1 


370 


2.1 


2012-06-03 


04:54 


F 


17.9 


0.356 


1 


370 


2.2 



-169°, 19° 
-40°, 8° 
169°, 25° 
172°, 25° 
158°, 27° 
109°, 28° 

-102°, 28° 
-59°, 28° 

-137°, 28° 
159°, 28° 

-19°, 10° 
25°, -71° 
-61°, -77° 
-57°, -78° 
-34°, -9° 
-83°, -9° 

65°, -9° 
107°, -8° 

29°, -8° 
-35°, -8° 



Table 3: Sub-observer points for all published visible- 
wavelength spectra of 1999 JU3. 

Columns in this table are the UT midpoint of observa- 



tion, reference [BOl = Binzel et al. (2001a); V08 = Vi- 



las ( 2008[ ); letters A-F indicate the corresponding spectra 



in Figure [T], solar phase angle (a), and the sub-observer 
longitude and latitude. These data are presented for both 
shape models of 1999 JU3. 
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Fig. 1. Normalized visible- wavelength reflectance spectra of 1999 JU3. The data are 
vertically offset for clarity. The top panel shows the three previous visible spectra 
with the month/year of observation and respective references BOl = Binzel et al. 



(2001a) and V08 = Vilas (2008). The July 2007 data is shown in grey to distinguish 



it from the adjacent spectrum. Our spectra are in the bottom panel and are chrono- 
logically ordered from bottom to top with spectra A and B taken on June 1, C on 
June 2, and D-F taken on June 3. The midpoint times for each of these observations 
are given in Table [3j The progressive decrease in signal-to-noise in spectra A-F is 
due to the opposition effect, increasing geocentric and heliocentric distances, and 
increasing background from a waxing moon. 
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Y-86029 



0.6 



0.5 



1.0 1.5 

Wavelength (microns) 
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Fig. 2. Composite visible (LDSS3) and near-infrared (SpeX and FIRE) spectrum 
of 1999 JU3. The grey region in the background represents the envelope for C-type 
spectra in the DeMeo et al. (2009) taxonomic system. In spite of the low signal- 
to-noise in the near-infrared and the presence of residual telluric features around 
1.4 and 1.9 fim, this spectrum shows no evidence for prominent absorption bands, 
though subtle features may be present around 0.9 fiui and short-ward of 0.7 //m. 
The two best spectral matches from our search through RELAB are a thermally 
altered sample of Murchison and the unusual CI chondrite Yamato 86029. These 
meteorite spectra have been normalized at 1 fj,m and then offset by -0.15 and -0.3 
units respectively. 
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Fig. 3. Rotational light curve of 1999 JU3. Filled circles represent the IMACS Mag- 
ellan data. Open diamonds represent our additional observations (Table [2]) with 
magnitudes scaled to the IMACS data. The solid line represents the predicted light 



curve from the 


Miiller et al. 


(2011) 


Kawakami et al. 


(2010 


). The data 



model, the dashed is the model prediction from 
have been phase folded to a synodic period of 



7.631 hours. A rotation phase of zero corresponds to the start of IMACS obser- 
vations on UT 2012-04-05 at 06:40 (JD 2456022.778). The phases accessed by the 
visible- wavelength spectroscopic observations are indicated by the letters A-F. The 
total time for each spectrum (4 x 180 seconds plus readout and telescope offsets) 
is comparable to the width of the letters. Spectroscopic coverage was obtained over 
approximately 60% of rotation phases at a longitudinal resolution of ~ 45°. 
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0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Wavelength (microns) 

Fig. 4. Visible-wavelength reflectance ratios of 1993 JUS relative to a median com- 
bination of all spectra. These reflectance ratios have been offset by mnltiples of 0.4 
units. Deviations from a slope-zero line (black curves) would indicate spectroscopic 
heterogeneity. No variability is detected within the S/N of the data, though slight 
(< 4%) deviations from a fiat line are seen at the shortest wavelengths in spectra 
B and C. 
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Fig. 5. Cylindrical projections of 1999 JU3's surface with circles indicating the sub- 



-observer latitudes and longitudes from the two available shape models ( Kawakami 



et al. , 2010 Miiller et al. 2011 ). The circles are not meant to indicate the physical ex- 
tent of the surface accessed by the observations. Black circles represent sub-observer 
points from our new 2012 observations and are labeled A-F. Grey circles represent 
sub-observer points for the observations from 1999 and 2007. The two spectra that 



display absorption features are indicated with hatched regions and labeled B ( Binzel 



et al. 


2001a) and V ( 


Vilas 


2008 



must be limited in extent to avoid overlapping those corresponding to featureless 
spectra. The loss of rotational information across epochs (see text) suggests that the 
locations of these regions relative to the 2012 points could lie anywhere along the 
horizontal bands. This demands that any heterogeneity on the surface be confined 
to regions no more than a few hundred meters in extent. 
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Fig. 6. SUPPLEMENTARY ONLINE FIGURE: Cylindrical projection of 1999 
JU3's surface with sub-observer hemispherical footprints indicating the regions ac- 
cessed by our spectra A-F. Projections are shown for both the Miiller et al. (2011) 
and Kawakami et al. (2010) shape models. The hemispherical regions are color 
coded for each spectrum as indicated by the letters in the top panel. The sub-ob- 
server hemispheres for spectra A and F are so similar that they have been combined 
into a single dark blue region. The sub-observer points in Figure 5 would lie at the 
center of their respective regions. The Vilas (2008) and Binzel et al. (2001a) hemi- 
spheres are not plotted, but it is clear that they would almost entirely overlap with 
the regions here. Our spectra obtained hemispherical averages over a majority of the 
asteroid's surface, independent of which shape model is assumed, thus leaving little 
room for heterogenous regions that could produce prominent spectral absorption 
bands. 
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